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ABSTRACT 

Galaxies are not distributed randomly in the cosmic web but are instead arranged in 
filaments and sheets surrounding cosmic voids. Observationally there is still no con- 
vincing evidence of a link between the properties of galaxies and their host structures. 
However, by the tidal torque theory (our understanding of the origin of galaxy angular 
momentum), such a link should exist. Using the presently largest spectroscopic galaxy 
redshift survey (SDSS) we study the connection between the spin axes of galaxies and 
the orientation of their host filaments. 

We use a three dimensional field of orientations to describe cosmic filaments. To 
restore the inclination angles of galaxies, we use a 3D photometric model of galaxies 
that gives these angles more accurately than traditional 2D models. 

We found evidence that the spin axes of bright spiral galaxies have a weak tendency 
to be aligned parallel to filaments. For elliptical/SO galaxies, we have a statistically 
significant result that the spin axes of ellipticals are aligned preferentially perpendic- 
ular to the host filaments; we show that this signal practically does not depend on the 
accuracy of the estimated inclination angles for elliptical galaxies. 

Key words: methods: statistical - galaxies: general - galaxies: statistics - galaxies: 
evolution - large-scale structure of Universe. 



1 INTRODUCTION 

Disc-dominated galaxies are the morphologi cally dominant 
class of galaxies in the present-day Universe (|Bamford et al.l 
2009). As these galaxies are rotationally supported, it is of 
vital importance to understand how disc galaxies acquire 
their angular momentum. The tidal-torque theory predicts 
alignment effects of disc galaxies, since the acquisition of the 
angular momentum is partially governed by environmental 
effects, such as tidal shearing produced by the neighbouring 
primordial matter distribution and the mo ment of inerti a 
of the forming protogalaxy (for a review see Schafcr 2009). 
Hence, testing the alignment of the orientation of galaxies 
with that of their surrounding environment (e.g., filaments) 
provides vital information about how galaxies form in the 
cosmic web. 

It is known that a ngular mom e nta o f nei g hbour 
galaxies are correlate d (|Slosar et all 120091 : iLeel l201ll : 
lAndrae fc Jahnkel201ll ), indicating that the environment in- 
fluences the acquisition of the angular momentum. However, 
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lAndrae fc Jahnkd i|201ll ) argue that this correlation is plau- 
sible but not statistically significant. 

During the last years, correlations between dark 
matter haloes and their host fil aments have be e n de- 
tected in Af-bod y simul ati ons (iHatton k. Nininl |200ll : 
[ Faltenbacher et all 120021: [Bailin fc Steinmet j 120051 : 



lAltav et all 120061: iBrunino et all 120071). More specifi - 



call y. lAragon-Calvo et all l|2007l ). iTrowland et~ai1 l|2012|} . 



and ICodis et al. I l|2012l ) showed that the orientation of the 
halo spin vector is mass-dependent. Spin axes of low-mass 
haloes tend to be aligned parallel to the filaments, whereas 
high mass haloes have an orthogonal alignment. In this 
picture, low-mass haloes form through the winding of flows 
in cosmic sheets/walls; hence at the intersection of these 
walls (in filaments), the spin axes of the haloes are parallel 
to the filaments. On the other hand, massive haloes are 
typically the product of mergers along such a filamentary 
structure, and thus their spin axes are perpendicular to the 
spine of the filaments. 

Previous studies based on observed galaxy catalogues 
have also indicated that spin axes of galaxies tend to 
be correlated with the structures in which they are em- 
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bedded dKashikawa fc Okamural Il992l ; iNavarro et all |2004 
iTruiillo et al.ll2006l ). These studies concentrate on sheets and 
voids and show that the spin axes of galaxies are preferen- 
tially oriented per pendicular t o the host structures. In a 
more recent study. Hones et all |2010l ) used a small sample 
of edge-on galaxies from the SDSS and showed that the spin 
axes of these galaxies are aligned perpendicular to the spine 
of the parent filament. These objects are also less bright, 
and so their results are in contradiction with the results of 
iV-body simulations. 

The picture is even more complicated since the spin 
axes of galaxies and their host h aloes are not strictly par- 
allel; Ivan den Bosch et al.l (|2002l ) find a median misalign- 
ment of approximately 30°. Furthermore, minor mergers 
can significantly di sturb the angular momenta of galaxies 
(|Moster et al.ll2010l ). 

To add even m ore controversy to the picture, 
ISlosar fc White! (|2009l ) claimed recently that in contrast to 
previous studies, there are no correlations between galax- 
ies and the v oids i n which they are located. Furthermore, 
IVarela et all (|2012h used the SDSS data and the morpho- 
logical classification from the Galaxy Zoo project and found 
that galaxy spins tend to be perpendicular to the void walls 
th ey are located in, whi ch is in contradiction with the results 
of lNavarro et al.l (120041 1 and lTruiillo et all (120061 ). Thus, the 
results remain contradictory, mostly because of the difficulty 
to measure the inclination angles of galaxies and/or to prop- 
erly define the large-scale filamentary structures. 

In the present study, we use the full SDSS spectro- 
scopic sample to study the correlation between the spin axes 
of galaxies and the orientation of their host filaments. We 
build a field of orientations for the filaments. We use a three- 
dimensional galaxy model to estimate the inclination angles 
of galaxies. To have good estimates of the inclination angles, 
we limit our sample strictly to galaxies with the bulge-to-disc 
ratio less than three: e.g. spirals and elliptical/SO galaxies. 

The structure of the present paper is as follows: in 
Sect. [5] we briefly describe the used data. In Sect. Owe de- 
scribe the three-dimensional galaxy model and explain how 
we selected the galaxy sample for the correlation study. The 
algorithm to build the filament orientation field is described 
in Sect. [4] and the correlation estimator is defined in Sect. [5] 
Finally, in Sect. [S] we give the results and discuss them. 

Throughout this paper we assume the following cos- 
mology: the Hubble constant Ho = 100/ikms _1 Mpc _1 , the 
matter de nsity f? m = 0.27 and the dark energy density 
n A = 0.73 (|Komatsu et al.ll201ll ). 



2 SDSS DATA 

Our present study is based on the SDSS DR8 ( Aihar a~et al.l 
l201ll ). We use only the main contiguous area of the sur- 
vey (the Legacy Sur vey) and the s p ectros copic galaxy sam- 
ple as compiled in iTempel et al.l |2012i ). The lower Pet- 
rosian magnitude limit for this sample is set to m r = 17.77, 
since for fainter galaxies, the spectroscopic sample is incom- 
plete l|Strauss et al.ll2002l ). To exclude the Local Superclus- 
ter from the sample, the lower CMB-corrected distance limit 
z = 0.009 was used. The upper limit was set to z = 0.2, 
since at larger distances the sample becomes very diluted. 
The sample includes 576493 galaxies. 



The DR8 sample is greatly improved over the previous 
data releases, in the sense that all the imaging data have 
been reprocessed using an improved sky-subtraction algo- 
rithm, a self-consistent photometric recalibration, and the 
flat-field determination. This improvement makes the DR8 
a good sample for detailed photometric modelling. 

For photometric modelling we used the atlas (pre- 
fix 'fpAtlas'), mask (prefix 'fM') and point-spread- function 
(prefix 'psField') images of the sample galaxies. The im- 
ages were downloaded from the SDSS Data Archive Server 
(DAS), and the SDSS software utilities readAtlasImages 
and readPSF were used for the atlas and PSF images, re- 
spectively. Masks were applied to the atlas images, using 
galaxy positions (given by rowc and cole in SDSS CAS). 

For photometric modelling, the initial guess for the cen- 
tre, position, and size of a galaxy is taken from the SDSS 
Catalog Archive Server (CAS), where each galaxy has been 
approximated by simple de Vaucouleurs and exponential 
models. To model the structural parameters, only the gri 
filters were used, since the uncertainties for the u and z 
imaging are the largest. 

To find the filamentary structure in the SDSS data, we 
have to suppress first the finger-of-god redshift distortions 
for groups. For th at we use the Frien ds-of- Friends (FoF) 
groups compiled in ITempel et all |2012l ): the details of the 
group finding algorithm are explained in iTago et al l (|2008l . 
l20ld ). We spherisize the groups using the rms sizes of galaxy 
gr oups in the sky and the ir rm s radial velocit i es as d escribed 
in lLiivamagi et al.l l|2012l) and ITempel et alj (|2012l ). Such a 
compression will lead to a better estimate of the density field 
and can help to find the real filamentary structure. 

To define the filamentary structure, we use a Cartesian 
grid based on the SDSS angular coordinates r\ and A, al- 
lowing the most efficient placing of the galaxy sample cone 
inside a box. The galaxy coordinates are calculated as fol- 
lows: 

x — — dgai sin A 

y = dgai cos A cos r\ (1) 
z = dg a i cos A sin 77, 

where d ga i is the finger-of-god suppressed co-moving dis- 
tance to a galaxy. 

The main purpose of this paper is to study the correla- 
tion between the spin vectors of galaxies and the spine/axis 
of filaments. Since the full angular momentum of galaxies 
cannot be estimated from photometric observations, we use 
the normalised spin vectors that describe the rotation axis 
of galaxy. Moreover, the spin vectors of galaxies (specifically, 
the inclination angles of galaxies) can in principle estimated 
only for disc-dominated galaxies: for pure elliptical galaxies 
the inclination angle cannot be found uniquely from pho- 
tometric observations. To select spiral galaxies and lentic- 
ular (SO) gala xies, we use the mo rphological classification 
as described in ITempel et al.l (|201ll ) , where we divide galax- 
ies into three classes: spirals, ellipticals, and with uncertain 
classification. This classification takes into account the SDSS 
model fits, apparent elliptic ities (and app a rent s izes), and 
different galaxy colours. In ITempel et al.l (|2012h we com- 
pared this cl assification with the cla ssification by the Galaxy 
Zoo project l|Huertas-Companv et al]|201ll ) and showed that 
they are in very good agreement. 



© 2012 RAS, MNRAS OOO.rflfm 



Spin alignment of galaxies in filaments 3 



So, for finding the filaments and their orientation we use 
the full SDSS spectroscopic sample. For correlation studies 
we restrict our sample to the galaxies with the bulge-to-disc 
ratio less than three (see Sect. 13.3)) . 



3 THREE-DIMENSIONAL PHOTOMETRIC 
MODELLING 

In this section we describe the 3D photometric model and 
describe in detail the sample of galaxies that we will use 
to study the correlations between galaxies and filaments. 
Full results of the photometric modelling of the SDSS main 
galaxy sample together with a catalogue will be published 
in a separate paper. 



3.1 Description of the three-dimensional 
photometric model of galaxies 

We use a three-dimensional model to describe the structure 
of a galaxy. The model is fully descri bed in ITempel et all 
i|2010l . l201ll ) (in the case of M31). In lTempel et all (|2012T ) 
we applied the model to a test sample of SDSS galaxies 
and showed that it can be used to automatically model the 
SDSS galaxies with a sufficiently good accuracy. Especially 
important for the present study is that the inclination angles 
of galaxies can be restored very accurately: for the majority 
of disc-dominated galaxies the differences are less than 5 
degrees; for galaxies with the bulge-to-disc ratio less than 
three the accuracy is better than 10-15 degrees. 

For consistency reasons we give a short description of 
the model below. 

The model galaxy is given as a superposition of its indi- 
vidual stellar components. In the present study we use a two 
component model: a spheroid + a disc. Both components are 
approximated by an ellipsoid of revolution with a constant 
axial ratio q; its spatial density distribution follows Einasto's 
law 



1(a) = 1(0) exp 



(—) 

\kaoJ 



l/N 



(2) 



where 1(0) = hL/ (4nqag) is the central density and L is the 
luminosity of the component; a — r 2 + z 2 /q 2 , where r and 
z are the cylindrical coordinates. We use the harmonic mean 
radius ao to describe the extent of the component; h and 
k are normalising parameters d epending on the st ructural 
parameter N (see appendix B of lTenies et al.lll994T ). 

The density distributions of both components are pro- 
jected along the line-of-sight and their sum yields the surface 
brightness distribution of the model 



2 00 

odei(x,y) = 2j2^- f 

.7 = 1 3 , 



lj(a)ada 



1/2 ' 



(3) 



where A 2 = x 2 + y 2 /Q 2 is the major semi-axis of the 
equidensity ellipse of the projected light distribution and 
Qj are the apparent axial ratios of the components (Q 2 = 
cos 2 i + q 2 sin 2 i) ; i is the inclination angle between the spin 
vector of the galaxy and the line-of-sight. The summation is 
over the two visible components (the spheroid and the disc). 
Equation ((3| gives the surface brightness of the galaxy for 
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Figure 1. Distribution of the inclination angles of galaxies. The 
green solid line shows the distribution for all galaxies; the blue 
dashed line shows the distribution for spiral galaxies . Spiral galax- 
ies we re selected based on the classification given in lTempel et al.l 
J201 lh . The peaks in the distribution are the remnant of the mod- 
elling procedure (see text for more information). 



all lines-of-sight that were used to create the model image of 
a galaxy. Before comparing model images with observations 
we convolved the model images with the SDSS point spread 
function (PSF) for a given galaxy. 

The correctness of the model fit was estimated by using 
the x 2 value, defined by 



X 



iVdof ^ ^2 



[Jobs (^,2/) ~ /model (Z, V)] 



a»(x,y) 2 



(4) 



where iVdof is the number of degrees of freedom in the fit; 
f^ bs (x,y) and fmodei( x iV) are the observed and modelled 
fluxes at the given pixel (x,y) and index v indica tes the 
filter (gri); <r(x, y) is the Poisson error at each pixel l|Howelll 
2006). The summation is taken over all filters (v) and all 
pixels of each galaxy as defined by the corresponding mask. 

To minimise the \ 2 we have used the downhill simplex 
metho d of Nelder and M ead from the Numerical Recipe li- 
brary l|Press et al.lll9 92). This method is quite simple and 
efficient in s earching large param eter spaces. The test sam- 
ple study in ITempel et al.l (|2012T ) showed that the method 
performs well and gives sufficiently acc urate results with a 
reasonable time. In lTempel et al.l l(2012h we also gave an ex- 
ample how this model performs for the real SDSS data; there 
are no systematic trends and the model performs well. 



3.2 Estimation of the inclination angles of galaxies 

Since our method to model galaxies depends slightly on the 
initial parameters, due to the degeneracy between various 
parameters, we cannot find all the parameters directly. The 
degeneracy is strongest for the inclination angles of galaxies, 
which are strongly degenerate with the thickness/ellipticity 
of the galaxy. 

To measure the inclination angle of a galaxy, we run the 
modelling code with ten different initial guesses for the incli- 
nation angle between 0° (face-on galaxies) and 90° (edge-on 
galaxies). This modelling was done for the gri filters sep- 
arately, in order to have three different guesses for the in- 
clination angle. Based on the modelling in the three filters, 
we choose the best inclination angle value. For the major- 
ity of disc-dominated galaxies these three filters give close 
inclinations as the best values. 

In the final modelling step, the gri filter data were used 
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Figure 2. The distribution of the bulge-to-disc ratios. The green 
solid line — all galaxies; the blue dashed line - spiral galaxies; the 
red dotted line - elliptical galaxies. The grey vertical dashed line 
shows the limit ibulge/^disc = 3 that we used to select galaxies 
for the correlation study. 



together with the best guess for the inclination angle. The 
final modelling gives us the structural parameters of galaxy 
components and the position and inclination angles of the 
galaxy - the orientation of its spin vector. The structural pa- 
rameters together with the morphological classification were 
used to select mainly disc-dominated or elliptical/SO galax- 
ies for the correlation study. 

Figure [T] shows the distribution of the inclination angles 
of galaxies for all galaxies and for spiral galaxies. Compared 
with the uniform distribution the number of both face-on 
and edge-on galaxies is smaller. This is caused by two rea- 
sons: firstly, the inclination angle is degenerate with the 
thickness of the disc, so it is hard to find exactly edge-on 
galaxies automatically; secondly, real galaxies are triaxial 
ellipsoids and even for exactly face-on galaxies, the model 
sees it as a slightly inclined galaxy. In the distribution, we 
also see the peaks of the initial inclination angle values (ten 
values from 0° to 90°). However, the deviation from the 
random distribution is not severe and it does not bias our 
correlation studies, it just adds some noise. 

If we look at the distribution of the inclination angles of 
spiral galaxies (Fig. [1} , we see that these galaxies are mostly 
edge-on. This is purely a selection effect: spiral galaxies are 
much easier to classify if they are edge-on. In our spiral 
galaxy sample, we miss many face-on spiral galaxies. In the 
correlation study we will take this selection effect into ac- 
count and it does not influence our results. 



3.3 The sample of spiral and elliptical/SO galaxies 

To study the correlation between galaxies and filaments we 
use spiral and elliptical/SO galaxies. The reason is that the 
inclination angles of galaxies can be adequately restored 
only for spiral/SO galaxies; for purely elliptical galaxies the 
inclination angles cannot be uniquely recovered. Our mor- 
phological classifica tion of galaxies was ini tially done in 
ITempel et ail l|201lf) . In ITempel et al.l l|2012h we combined 
this classification with the morphological classificatio n by 
the Galaxy Zoo project (|Huertas-Companv et al]|201lf) and 
improved our classification. 

In this paper, we use the results of our galaxy modelling 
and restrict our galaxy sample to the objects with the bulge- 
to-disc ratio less than three. Figure [2] shows the distribution 
of the bulge-to-disc ratios for all, spiral, and elliptical galax- 
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Figure 3. Two dimensional projection of a thin cylinder with its 
shadow within a pattern of galaxies. 



ies. Most spiral galaxies have the bulge-to-disc ratio less than 
two, and most ellipticals have the bulge-to-disc ratio greater 
than one. 

To study the orientation correlation between galaxies 
and filaments we use two samples. Firstly, we use a spiral 
galaxy sample with Lbuige/idisc < 2. The second sample is 
an elliptical galaxy sample (by our morphological classifi- 
cation): for this sample we also use only the galaxies with 
ibuige/idiac < 3, representing mostly lenticular galaxies. 

As a final step, we eliminate from our sample these 
galaxies which are close to the borders of the survey since 
we cannot find reliable filament orientation there. Since the 
minimum filament length in our study is 5/i _1 Mpc (see 
Sect. we set the limiting distance from the survey border 
to half of it (2.5 /t~ 1 Mpc). We use the spatial sample mask 
as defined in Martmez et ail (|2009T ). 

As a result, our final sample that we use for the corre- 
lation study includes 231544 spiral galaxies and 26140 ellip- 
tical/SO galaxies. 



4 FINDING THE FILAMENTARY 
STRUCTURE 

4.1 Detecting filaments and their orientation 

In this section we describe the main tools we use to study 
the large-scale filaments orientations. Our previous papers 
i|Stoica et al.l 120071 . |2010D concentrate on detecting a gen- 
eral filamentary pattern. In this paper, our main motivation 
is to detect the filament orientation. To do that, we mod- 
ify/redefine the underlying model and use this model to de- 
tect the orientation of fila ments at galaxy locations. This 
approach is close to that of Irlill et al.l (|201ll ) who show how 
to estimate the field of orientations from the point pattern 
data. 



4-1.1 Elementary cylinder 

We model the filamentary structure by a network of thin 
connected cylinders. We assume that locally, galaxies may 
be grouped together inside a rather small thin cylinder. We 
also assume that such small cylinders may combine to extend 
a filament if neighbouring cylinders are aligned in similar di- 
rections. We used suc h a cylinder netwo r k to build filaments 
in our earlier papers (|Stoica et al 1 l2007l . l2010T l. In these pa- 
pers the cylinders were used to build an object point process 
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Figure 4. A sample of filaments and their orientations. Three various regions from the SDSS volume (units are in /i -1 Mpc) are shown. 
Blue dots show the location of galaxies. Grey lines show the location of probed cylinders (see text for more information). Red lines 
around galaxies show the filaments and their orientation. The figures are rotated so that the longest filaments are well seen. 



that delineated filaments; here we filter the point process 
given by galaxy positions by the cylindrical kernel to esti- 
mate the orientation field. This i s similar to the work d one in 
Aragon-Calyo et al.l (|2007l , l201Ch . ISousbie et ail j2008h . and 
Smith et al.l ((2012J), where the starting point is to determine 
a smoothed density field and to use its structure to define 
filaments. The main difference with our approach is that the 
kernel used to detect filaments in these papers is, by neces- 
sity, isotropic, while our kernel is anisotropic, cylindrical, 
and follows the local anisotropy of the galaxy distribution. 

A cylinder y is an object characterised by its centre k 
and its shape parameters. The shape parameters of a cylin- 
der are the radius r, the height h, and the orientation vec- 
tor uj. We consider the radius and the height parameters as 
fixed, selecting r — 0.5 h~ 1 Mpc, and h — 5.0 h~ Mpc. These 
values fit well the observed filamentary structure. The ori- 
entation vector parameters uj = <f>(ri, r) are uniformly dis- 
tributed on M = [0, 27t] x [0, 1], such that 

uj — (^\J I — t 2 cos(r/), \/\ — t 2 sin(r;), . (5) 

We denote the shape of a cylinder by s(y) = s(k, r, h, uj). 

For every elementary cylinder y, we introduce a bigger 
cylinder with exactly the same parameters as the elementary 
cylinder, except for the radius, which equals 2r. Let s(y) 
be the shadow of s(y) obtained by the subtraction of the 
initial cylinder from the bigger cylinder. The cylinder and 
its shadow are shown in Fig. [3] Then, each cylinder is divided 
in three equal volumes along its main symmetry axis, and 
we denote by si(y), S2(y), and ss(y) their corresponding 
shapes, and by §i(y), S2(y), and ss(y) their corresponding 
shadow shapes. 

We describe the galaxy data as the point field given by 
the set of the galaxy positions, d = {di, di, ■ ■ ■ , d n }. To find 
a cylinder that is allowed by the data, we check a number 
of conditions. We have found that it is advantageous to use 
slightly more strict criter i a tha n we used in our previous 
papers l|Stoica et al.ll2007l . [2oTol '). 

First, we demand that the density of galaxies inside s(y) 
is sufficiently higher than the density of galaxies in s(y), 
formalised as follows: 

n(d n s(y))/v(s(y)) > 2 ■ n(d n s(y))/K%))> (6) 



where n(dn s(y)) and n(dn s(y)) axe the numbers of galax- 
ies covered by the cylinder and its shadow, and v(s(y)) and 
v(s(y)) are the volumes of the cylinder and its shadow, re- 
spectively. So, we require that the density in the cylinder is 
at least two times larger than in its shadow. Additionally, 
for every cylinder element Si(y) we require that the density 
in the element should be larger than in its shadow: 

n{d n Si(v))/y( Si (y)) > n(d n h{y)) /v{Hv)), (7) 

where i £ (1, 2, 3) denotes the cylinder parts. 

To ensure that the cylinder is approximately uniformly 
populated, we require 

max(n(dn s 1 (y),n(dn s 2 (y),n(dn s 3 (y)) < g 
min(n(d n si(y), n(d n S2(y), n(d n S3(j/)) ~ 

This criteria also helps to forbid cylinders partially inside 
galaxy groups. 

If all those criteria are satisfied, then we accept the 
cylinder. These criteria are empirical priors, chosen by obser- 
vational experience. The most common observed filaments 
are filaments that connect groups, and the thickness of such 
filaments is about 1 /i _1 Mpc. This is the reason, why we 
chose r = 0.5 /i _1 Mpc: such filaments should be the sites of 
galaxy formation. 

If the above rules are satisfied, we calculate the weighted 
number of galaxies inside a cylinder (its richness): 

n 

where r denotes the cylinder radius, Di denotes the distance 
of a galaxy from the axis of the cylinder, and summation is 
over all galaxies inside the cylinder (n = n(d n s(y))). If 
the rules listed above are not satisfied, we set the richness 
(the intensity of the orientation field) at the centre of the 
cylinder to zero. 

The richness R(s(y)) is used to estimate the goodness 
of a cylinder: it is determined by the number density inside 
a cylinder, and weighting helps to find the best orientation 
for a cylinder. 
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4-1.2 Finding the filaments and their orientation 

To study the correlation between galaxies and filaments, we 
need to know the filament orientation at the location of a 
galaxy (if the galaxy is in a filament). To define the filaments 
and their orientations, we build a filament orientation field 
that has direction an d intensity at ev ery point of the sam- 
ple volume (see, e.g., iHill et al] l201ll ). For that we do the 
following. 

(i) We use a local morphological filter to detect the direc- 
tion and an intensity of the filament orientation field. This 
filter is based on the elementary cylinder defined in the pre- 
vious Section. This cylinder is defined only in the regions 
(in location and orientation space) that satisfy the criteria 
given above and this defines the mask for the field; out- 
side the mask, the intensity of the field is zero. The richness 
R(s(y)) defines the intensity of the field. 

(ii) In order to apply this filter, we discretise the loca- 
tion space, covering it with cells of size of l/i _1 Mpc. Our 
test calculations showed that this spacing is small enough to 
sample the direction field. We assume that only one cylin- 
der (centre of the cylinder) can be located in such a cell, but 
clearly a cell can be touched by several cylinders. 

(iii) To estimate the orientation field, we maximise the 
richness R(s(y)) in every cell inside the mask volume; we 
accept only those cylinders, which satisfy the criteria given 
in Sect. BXTl 

(iv) The maximisation of R(s(y)) is done using a fol- 
lowing maximisation procedure: inside each location cell, 
the centre position of the cylinder together with its corre- 
sponding orientation parameters are moved stochastically 
till its richness is maximised. For that we use multidi- 
mensional nested sampling, applying t he Multinest tool 
ijFeroz fc Hobsonll2008l : lFeroz et al'jl20Qg| ). This is a parallel 
algorithm and maximises the richness R(s(y)) in reasonable 
time. 

Basically, this is cylinder based data filtering and those 
accepted cylinders define an intensity map for filaments; the 
cylinder orientations define an orientation map for filaments. 
Visually, the cylinders are located in beams along a filament. 
These maps allow to find the number of the cylinders touch- 
ing a galaxy and this allows to decide whether a galaxy 
belongs to the filamentary structure or rather to a cluster 
(crossing points of filaments). 

For every galaxy, in order to find the filament orien- 
tation in its location, we firstly find all the cylinders that 
touch (include) the galaxy. If there are no such cylinders, 
the galaxy is not located in a filament. If at least 2/3 of the 
touched cylinders have roughly the same orientation (the an- 
gle with the mean direction of all cylinders is less than 20°), 
we say that there is a filament: to be able to apply the 2/3 
rule, there must be at least three cylinders that touch the 
galaxy. If there is no such prominent beam, then probably 
it is just a crossing point of filaments, where the orientation 
is unclear. We reject those galaxies and do not use them for 
our correlation studies. 

If there is a filament, we use those cylinders that helped 
us to decide that there is a filament (that have the angles 
with the mean direction less than 20°) and average their 
directions to define the orientation of their beam (the ori- 
entation of the filament at that point). We use the cylinder 



All galaxies 
Spiral galaxies 
Elliptical galaxies 



50 100 150 200 250 300 350 400 450 500 
Distance [Mpc h 1 ] 

Figure 5. The fraction of galaxies in filaments for all, spiral, and 
elliptical galaxies as a function of distance from us. 




0.4 0.6 

cos(9) (line-of-sight) 

Figure 6. The distribution of the angle between the filaments 
and the line-of-sight: cos 9 = 1 for the filaments parallel with 
the line-of-sight; cos 6 = for the filaments perpendicular to the 
line-of-sight. 



richness R(s(y)) as the weight for averaging. This algorithm 
gives sensible filament orientations. We use only those galax- 
ies that belong to the filamentary structure to study the 
correlation between the orientation of galaxies and of the 
filament orientation field. 

The last operations, leading to the reliable filament ori- 
entation value at a position o f a galaxy, is an example of a 
morphological filter (see, e.g.. ISerralll982r i. 



4.2 An example of filaments 

Figure [4] shows an example of filaments in three different 
parts of the SDSS volume. The cubes are rotated so that the 
best filaments are well seen. We see that the algorithm finds 
good filaments and the filament orientation is also plausible. 
A closer look reveals that filaments intersect in the location 
of groups. 

Grey lines in Fig.Ushow the allowed cylinders and dark 
red lines around galaxies mark the galaxies which are in fil- 
aments, together with filament orientation. Initially there 
are a few cylinders that do not belong to filaments (that 
are single and that are discarded when calculating filament 
orientation), but if we test the galaxies against the cylinder 
pattern, we can eliminate those cylinders (it is well seen in 
upper right part in central panel of Fig. [4|. Our method to 
find the filaments depends on the number density of objects: 
if the number density is too low, many of the filaments are 
undetectable. Figure [5] shows the fraction of galaxies in fila- 
ments as a function of distance from the observer. It is seen 
that up to 200 /i _1 Mpc, the fraction is almost constant and 
decreases afterwards. The fraction of galaxies in filaments 
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is almost the same for spiral and elliptical galaxies, being 
slightly larger for ellipticals. 

Figure [6] shows the distribution of the angles between 
the filaments and the lines-of-sight. It is seen that our fila- 
ment finder finds less filaments parallel to the line-of-sight. 
The reason for this is unclear, but there are two reasonable 
possibilities. Firstly, since we suppress the finger-of-god red- 
shift distortions, we may over-suppress the groups; or we 
may include galaxies in the front or back of the groups that 
actually belong to filaments, into groups. The second reason 
lies in large-scale cosmological redshift distortions that com- 
press sup erclusters along th e line-of-sight in velocity space 
(see, e.g-. lPraton et al.lll997T l. For correlation study, we take 
this selection effect into account. 



5 ESTIMATING THE GALAXY-FILAMENT 
CORRELATION 

To quantify the impact of the filament environment on 
galaxies, we start with our estimate of the filament orien- 
tation field, found using all galaxies in the sample. Using 
this, we study how the spin vectors of different galaxy sub- 
samples are correlated with the orientation of the filament 
in which the galaxies reside. 

5.1 Computation of the alignment signal 

To estimate the orientation of galaxies relative to filaments, 
we compute the probability distribution function P(cosO) 
(alignment signal), where 8 is the angle between the rota- 
tion axis of a galaxy and the spine/orientation of the fil- 
ament. The quantity cos 8 is obtained as a scalar product 
between the unit vectors of spin axes of galaxies (s ga i axy ) 
and orientation vectors of filaments (r n i amont ) 

COS 8 = | S ga l axy ■ r filament | • (10) 

Here, cos 8 is restricted to the range [0,1], and cos# = 1 
implies that the rotation axis of galaxy is parallel to the 
filament while cos 8 = indicates the perpendicular orien- 
tation. 

The probability distribution function should be com- 
pared with the null hypothesis of random mutual orienta- 
tion of galaxies and filaments. We have seen that neither 
our galaxies nor filaments have random orientations due to 
selection effects (see Figs. [T]and[|)]). We use Monte-Carlo ap- 
proximation to estimate the orientation distribution for the 
case where there are no correlations between the orientations 
of filaments and galaxies, a nd to find the confi dence intervals 
for this estimate (see, e.g.. Illlian et al.ll2008l ). For that, we 
generate 1000 randomised samples in which the orientations 
(spin vectors) of galaxies are kept fixed, but galaxy locations 
are randomly changed between each other. This gives a new 
filament orientation (in respect to the line-of-sight) at the 
location of a galaxy, and a new mutual orientation angle. 
For each of our subsamples we find the median of the align- 
ment distribution for the null hypothesis together with its 
95% confidence limits. This distribution takes into account 
both the observational selection effects and the influence of 
the sample size on the estimate of the distribution. These 
data are shown in every subsequent figure (Figs. [7TTT01). for 
comparison with the observational results. 



5.2 Spin vectors of galaxies 

To compute the angle 8 we need to define the spin vector of 
a galaxy (s ga iaxy) and the orientation vector for a filament 
(rfliament)- The last one is uniquely defined as described in 
Sect. [4] For galaxies, we need to know the position angle 
(orientation in the sky) of a galaxy and its inclination angle. 
Both of them are estimated by galaxy modelling as described 
in Sect. O but since we do not know which side of the galaxy 
is closer to us, we have two spin vectors for each galaxy 
(actually four, but we are neglecting the direction of the 
spin). 

Different approaches have been used to deal with the in- 
determinacy of the spin vectors of galaxies. Since for edge-on 
and face-on galaxies, the indeterminacy is i rrelevant, sev- 
eral s tudies have used onl y those galaxies l|Truiillo et al.l 
120061 ; ISlosar fc White] |2009| ). The disadvantage of this ap- 
proach is the fact that the number of galaxies in the sample 
is greatly re duced. Another possi bility is to use just one 
spin vector dLee &: Erdogdul | 2007|) or both simu lt aneously 
|Kashikawa fc Okamuralll992l ; IVarela et al.ll2012l ). However, 
this approach dilutes the strength of the measured signal. 

Since we use the statistical approach to calculate the 
alignme nt signal, we decide d to use both spin vectors of 
galaxies. IVarela et al.l (|2012h also tested this approach with 
several Monte Carlo simulations and showed that the proce- 
dure recovers correctly the probability distribution function. 



6 RESULTS AND DISCUSSION 

In Table [T] we show the number of galaxies, the average of 
cos#, and Kolmogorov-Smirnov (KS) test probability that 
there is no correlation between galaxies and filaments for 
subsamples presented in this Section. For the KS test, we 
compared the observed distribution of mutual orientations 
with the mean of the randomised distributions of mutual 
orientations (P(cos #) rn d), as explained in Sect. I5TT1 In or- 
der to compare our results with previous work, we also cal- 
culated the average (cos#) of the orientation distribution. 
As the probability density for observed angles P(cos#) b s 
is conditional on the distribution of random orientations 

P(C0S 6>)rnd, 

P(cos 8) obs = /(cos 6»)P(cos 8) lnd , (11) 

we find the estimate /(cos 8) from the above formula and 
use it to estimate the mean (cos 8). 

6.1 Spiral galaxies 

Figure [7] shows the probability distribution P(cos8) for the 
mutual angle 8 between the spin axes of spiral galaxies and 
the orientation vectors of filaments. The grey shaded area 
corresponds to the random distribution (95% confidence re- 
gion) and the dark red line shows the real signal. It is seen 
that for all spiral galaxies, the real signal lies within the 95% 
confidence region for the random distribution. The KS test 
also confirms that the observed signal does not differ from 
the random distribution. 

The lower panel in Fig.[7|shows P(cos 8) for bright spiral 
galaxies (M r < —20.7 mag). For bright spiral galaxies, there 
are clearly more galaxies which have the spin axes parallel 
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Table 1. Number of galaxies, average cos 9, and the KS test 
probability that the sample is drawn from a uniform mutual ori- 
entation for our subsamples of galaxies. 
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refers to the spiral sample, where the inclination is calculated 
from Eq. JT2J. W refers to the elliptical/SO sample where the 
inclination angles are randomised (see below). refers to the 
elliptical/SO sample where the position angles are randomised (see 
below). 




cos(8) 

Figure 7. The orientation probability distribution for all (upper 
panel) and bright (lower panel) spiral galaxies. The black line and 
the grey filled region show the null hypothesis together with its 
95% confidence limits. The solid red line shows the true alignment 
distribution. 



to filaments. This is confirmed by the KS test and (cos 6*) 
(see Table [TJ). However, since we can model brighter galax- 
ies better (and can better estimate the inclination angles), 
we cannot definitely say that the correlation is intrinsically 

stronger for brighter gal axies. 

Many studies (e.g. ISlosar fc White] 120091 : IVarela et aD 
l20T2j) use the apparent axial ratio (b/a) and an assumed 
disc flatness (/) to restore the inclination angles i of galaxies. 
They use the formula 



(b/a) 2 - f 
1 _ p 



(12) 



For values of b/a < /, the angle is set to 90°. The flatness / 
depends on the morphological type of galaxies, but usually 
an average value is used. 

To test how this simplified equation can be used to 
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Figure 8. The orientation probability distribution for bright spi- 
ral galaxies, where the inclination angles are estimated using 
Eq. U2t . The black line and the grey filled region show the null 
hypothesis together with its 95% confidence limits. The solid red 
line shows the true alignment distribution. 
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Figure 9. The orientation probability distribution for all (upper 
panel) and bright (lower panel) elliptical/SO galaxies. The black 
line and the grey filled region show the null hypothesis together 
with its 95% confidence limit. The solid red line shows the true 
alignment distribution. 



study the correlation between the orientation of galaxies and 
filaments, we calculated the inclination angles u sing Eq . (112I) 
and a ssuming constant disc flatness / = 0.14 (|Varela et al.l 
2012), which is in good agreement with the average value 
from our 3D modelling. 

Figure [5] shows the correlation between the spin axes 
of spiral galaxies and filaments, when the inclination angles 
were estimated using Eq. (|12[) . It is seen that the correla- 
tion is weaker, showing that the simplified equation is not a 
good tool to estimate the inclination angles of galaxies. As 
a conclusion, detailed photometric modelling is necessary to 
restore the inclination angles accurately. 



6.2 Elliptical galaxies 

We also study the correlation for elliptical/SO galaxies. For 
that we restrict our sample to galaxies where the bulge-to- 
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Figure 10. The orientation probability distribution for bright 
elliptical galaxies. The upper panel shows the distribution when 
inclination angles are randomised. The lower panel is for ran- 
domised position angles. 



disc ratio is less than three, because for higher bulge-to-disc 
ratios, the inclination angle cannot be restored. 

Figure [9] shows the cos 9 probability distribution for all 
and for bright elliptical galaxies. As for spiral galaxies, the 
correlation signal is larger for brighter galaxies. (The KS 
test p-values are 3.6 • 10 -7 and 1.5 • 10 -9 , respectively.) 

Actually, if we use only faint elliptical galaxies, the sig- 
nal is almost lost. There might be two reasons: firstly, the 
correlation signal is weaker for fainter galaxies; secondly, this 
is purely a selection effect. For fainter galaxies, the mor- 
phological classification is less accurate and also the PSF 
makes galaxies rounder and the position/inclination angles 
are much harder to estimate. Compared with spiral galaxies, 
the correlation signal is much stronger. 

Since the inclination angles for elliptical galaxies cannot 
be restored with high accuracy, we studied how the corre- 
lation signal depends on inclination angles. For this tests, 
we used the bright elliptical/SO sample. In the first test, we 
randomised the inclination angles of galaxies. The results 
are shown in the upper panel of Fig. 1101 It is seen that 
for randomised inclination angles, the correlation signal re- 
mains, meaning that the signal comes from filaments which 
are perpendicular to the line-of-sight. If the spin axis of a 
galaxy is perpendicular to such a filament, changing the in- 
clination angle does not change the angle between the spin 
of the galaxy and the filament. 

In the second test, we kept the inclination angles fixed 
and randomised the position angles of galaxies. The results 
are shown in the lower panel of Fig. 1101 The correlation sig- 
nal is much weaker in this case, indicating that the inclina- 
tion angles are not accurately restored. However, the corre- 
lation signal is still visible (the KS test p- value is 9.1 • 10 -4 ), 
meaning that the modelled inclination angles are statisti- 
cally close to the real ones. 

As a conclusion, the correlation signals for elliptical 



galaxies can also be calculated, regardless of the inaccuracy 
in the determination of the inclination angle. Moreover, the 
correlation signal is so strong, that it is visible even if the 
inclination angles are randomised. 

6.3 Comparison with TV-body simulations and 
previous observational detections 

Many TV-body simulations indicate that the correla- 
tion of spin axes of dark mat ter haloes with filaments 
(sheets) is halo mass dependent llAragon-Calvo et al. 2007; 
IZhang et ad1l2009l ; iTrowland et al.ll2012l ; ICodis et altolj ): 
spin axes of low-mass haloes tend to lie parallel to the fil- 
aments while high-mass ha loes have an ortho gonal align- 
ment. To be more specific, lHahn et al.l (|2010h studied es- 
pecially the orientation of disc galaxies in a hydrodynamic 
AMR simulation. They showed that the most massive galaxy 
discs have spins preferentially aligned so as to point along 
their host filam e nts. T he same conclusion was reached in 
iLibeskind et all l|2012l) . who claimed that low-mass haloes 
(that mostly host spirals) tend to be aligned with filaments. 
Both papers are in good agreement with our study, where we 
found that especially for brighter disc galaxies, the galaxy 
spins tend to be aligned parallel to the filaments. 

To compare TV-body simulations with our findings, 
we link low-mass haloes with spiral galaxies and high- 
mass haloes with elliptical galaxies. The probable formation 
mechanisms for low- and high-mass haloes support these 
connections. The spin directions of low-mass haloes result 
from the winding flows from the sheets into filaments. Such 
a process naturally produces a net halo spin parallel to the 
filaments. Since those haloes are not influenced by many 
mergers, they are natural hosts for spiral galaxies. 

On the other hand, high-mass haloes are significantly 
affected by mergers occurring mostly along the filaments. 
Therefore they acquire a spin which is preferentially per- 
pendicular to filaments. These merger-produced haloes host 
also elliptical galaxies which form through many mergers. 

Contrary to TV-body simulations, where there are many 
papers devoted to the study of the correlation between the 
orientation of galaxies and their host structures, there is 
little observational evidence for that. The reasons for that 
are, however, clear: it is hard to estimate the spin axes of 
galaxies and to identify the structures where galaxies are 
located (due to observational effects). Regardless of these 
difficulties, the re are still some pa pers devoted to this prob- 
lem. Recently, Ijones et alj (|2010l ) selected 69 well-defined 
edge-on galaxies and used the Multi-scale Morphology Fil- 
ter (MMF) technique to detect filaments. They found that 
less intrinsically bright galaxies tend to lie perpendicular 
to hos t filaments. Based on the study of voids IVarela et al] 
1120121) found that spin axes of disc galaxies around the voids 
are preferentially directed toward the centres of the voids. 
Since filaments are located around the voids, the spin axes 
are also perpendicular to filaments. 

Both these findings contradict with our result s . How - 
ever, our prese n t wor k agrees with iTruiillo et al. I (|2006l ); 
iLee fc Erdogdul (|2007h . where the authors found that spin 
axes of spiral galaxies are preferentially in the plane of the 
cosmic voids and parallel to filaments. In short, the observa- 
tional findings seem to be in disagreement with each other 
and the results remain inconclusive. The reasons for that are 
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still unclear. The biggest problem with observational find- 
ings is the statistical significance of the results and until this 
improves, we cannot say anything convincing. 

Our results support the tidal torque scenario where the 
galaxy spins indeed originated from the initial tidal inter- 
action with the surrounding matter. Moreover, our results 
are supported by recent high-resolution i V-body simulations 
l|Trowland et al.ll2012l ; ICodis et alj|2012h . 



7 CONCLUDING REMARKS 

We have examined the orientation of spin axes of galaxies 
with respect to the cosmic filaments. We used 3D photomet- 
ric modelling of galaxies to restore the inclination angles of 
galaxies with a sufficiently high accuracy, especially for disc- 
dominated galaxies. The alignment between galaxy spins 
and the axis of filaments was characterised by the shape 
of the probability distribution of cos# where 6 is the an- 
gle between the two vectors. We studied the correlation for 
spirals and ellipticals separately. 

We found a strong and significant correlation between 
the spin axes of elliptical galaxies and filaments; these galax- 
ies tend to be aligned perpendicular to filaments, whereas 
more luminous ellipticals have a stronger orthogonal align- 
ment. We showed that this finding is not influenced the fact 
that the inclination angles for elliptical galaxies cannot be 
estimated accurately. 

Within the caveats of a limited statistical significance, 
we found that the spin axes of bright spiral galaxies tend 
to be aligned parallel to the host filaments. The significant 
alignment we have found is due to finding the inclination 
angles of galaxies by 3D modelling: using a simple inclination 
angle estimate by the visible axial ratio produces a weaker 
correlation. This emphasises the benefits of having a good 
estimate for inclination angles when studying the correlation 
of galaxies with the cosmic web. 

Our findings for spiral and elliptical galaxies are in 
agreement with the recent results of iV-body simulations, 
suggesting that spirals form through peaceful accretion of 
matter and ellipticals are the results of mergers that mostly 
occur along the filaments. 

We hope that our results provide a better understanding 
of the formation of galaxies in the context of their host large- 
scale structures, which is one of key questions for galaxy for- 
mation theory. A precise knowledge of correlated galaxy ori- 
entations is of high importance also for future high-precision 
weak lensing studies of dark energy. So, it is of high impor- 
tance to improve the filament finding algorithms and the 3D 
modelling of galaxies to have a better observational insight 
into this problem. Future work is planned in this direction. 
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